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Infinite-range spin-glass models with Levy-distributed interactions show a freezing transition sim- 
ilar to disordered spin systems on finite connectivity random graphs. It is shown that despite 
diverging moments of the local field distribution this transition can be analyzed within the replica 
approach by working at imaginary temperature and using a variant of the replica method developed 
for diluted systems and optimization problems. The replica-symmetric self-consistent equation for 
the distribution of local fields illustrates how the long tail in the distribution of coupling strengths 
gives rise to a significant fraction of strong bonds per spin which form a percolating backbone at 
the transition temperature. 
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I. INTRODUCTION 

Spin-glasses are model systems of statistical mechanics in which simple degrees of freedom interact via random 
couplings drawn from a given probability distribution [lj. The ensuing interplay between disorder and frustration 
gives rise to peculiar static and dynamic properties which made spin-glasses paradigms for complex systems with 
competing interactions. In this way the concepts and techniques developed for their theoretical understanding Q 
became invaluable tools in the quantitative analysis of problems originating from such diverse fields as algorithmic 
complexity 0, 0, [H] , game theory @, 0] , artificial neural networks Q , and cryptography . 

In the present note we study a spin-glass for which the couplings strengths are drawn from a Levy-distribution. 
The main characteristic of these distributions are power-law tails resulting in diverging moments. Compared to 
the extensively studied spin-glass models with Gaussian or other finite moment distributions [ll|, [H, EH Levy- 
distributed couplings are interesting for several reasons. On the one hand the comparatively large fraction of strong 
bonds gives rise to a mechanism for the glass transition which is different from the usual scenario. On the other hand 
these systems pose new challenges to the theoretical analysis because the diverging second moments invalidate the 
central limit theorem which is at the bottom of many mean- field techniques. Related issues of interest include the 
spectral theory of random matrices with Levy-distributed entries [l4l . [l5| and relaxation and transport on scale-free 
networks [16]. It is also interesting to note that the characteristic properties of the Cauchy-distribution have recently 
enabled progress in the mathematically rigorous analysis of matrix games with random pay-off matrices [l7j . 

The model considered below with the help of the replica method was analyzed previously by Cizeau and Bouchaud 
(CB) using the cavity approach jl8j |. In their paper CB remark that they resorted to the cavity method because 
they were not able to make progress within the replica framework. This might have been caused by the fact that 
at that time the central quantity in the replica method was the second moment of the local field distribution, the 
so-called Edwards- Anderson parameter (l9| . which for Levy-distributed couplings is likely to diverge. Later a variant 
of the replica method was developed to deal with non-Gaussian local field distributions characteristic for diluted spin 
glasses and complex optimization problems [20] . Until now this approach was used only in situations where the local 
field distribution is inadequately characterized by its second moment alone and higher moments of the distribution 
are needed for a complete description. Here we show that the method may also be adapted to situations where the 
moments not even exist. 

II. THE MODEL 

We consider a system of N Ising spins Si = ±1, i = 1,...,N with Hamiltonian 

H({S i }) = -j± /s y £j ij S i S j , (1) 
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where the sum is over all pairs of spins. The couplings Jp = Jji are i.i.d. random variables drawn from a Levy 
distribution P a (J) defined by its characteristic function [23] 

P a (k) := J dJ e- lkJ P a (J) = e^ |fc| ° (2) 

with the real parameter a € (0,2]. The thermodynamic properties of the system are described by the ensemble 
averaged free energy 



with the partition function 



/(/?):=- lira m ^Z(p) , (3) 



Z(f3) := eM-PH({Si})) • (4) 



Here (3 denotes the inverse temperature and the overbar stands for the average over the random couplings J, 
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III. REPLICA THEORY 

To calculate the average in ([3]) we employ the replica trick [l9| 

7n _ 2 

In Z = lim . (5) 

ra->0 n 

As usual we aim at calculating Z n for integer n by replicating the system n times, {Si} i— > a = 1, ...,n, and 

then try to continue the results to real n in order to perform the limit n — > 0. 

Due to the algebraic decay P a (J) ~ | J| _Q_1 of the distribution P a (J) for large |J| the average Z n ((3) does not exist 
for real (3. On the other hand, for a purely imaginary temperature, /3 — —ik, k £ R, we find from the very definition 
of P a (J), cf. @ 



zH-ik) = x: cxp ( - ^ 2 1 e + • (6) 
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Note that the replica Hamiltonian is extensive which justifies a-posteriori the scaling of the interaction strengths 
with TV used in (Q]). The determination of Z n can be reduced to an effective single site problem by introducing the 
distributions 

i 

where S — {S a } stands for a spin vector with n components. We find after standard manipulations (2p| | 

Z^k) = J\[dc\S)5{Y J c{S)-l)eKv(-N[Y J c{S)\nc{S) + ^ . (8) 

s s s s,s> 

In the thermodynamic limit, N — > oo, the integral in ((8|) can be calculated by the saddle-point method. The 
corresponding self-consistent equation for c(<?) has the form 

c(a) = A(») cxp ( - \k\ a C 0)\S- <t\ q ) , (9) 

s 

where the Lagrange parameter A(n) enforces the constraint ^2sc(S) — 1 resulting from |7[). 
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IV. REPLICA SYMMETRY 



Within the replica symmetric approximation one assumes that the solution of ([9]) is symmetric under permutations 
of the replica indices implying that the saddle-point value of c(S) depends only on the sum, s :— ^2 a S a , of the 
components of the vector S. It is then convenient to determine the distribution of local magnetic fields P(h) from its 
relation to c(s) as given by [20| 



:(«) = J dh P{h) e- lkhs P(h) = J 7^ e lsh c(£) 



(10) 



Note that the P(h) defined in this way is normalized only after the limit n — ► is taken. The distribution of local 
magnetic fields is equivalent to the free energy ,f{(3) since all thermodynamic properties may be derived from suitable 
averages with P(h) [2l| . 

To get an equation for P(h) from ([9]) we need to calculate 

^e- lkhs \S ■ S\ a = J ^^\r\ a e ir rJ2exp(-ikhs-irS-a) (11) 
s s 

<t^L | r |« e irf ^ J[ exp ( - iS a (kh + f <r Q )) (12) 

s a 

' lr ' lr \r\ a e lrf [2cos{kh + f)]^ [2 cos(fc/i - r )] ^ (13) 



2vr 

cos(kh + f) 



drdf 
Zir 



cos(kh — f) 



(14) 



where the limit n — > was performed in the last line and a := a a . Using A(n) — > 1 for n — > |20| we therefore 
find from ([5]) in the replica symmetric approximation 

Using this result in (|10[) and performing the transformations r i— > r/fc, f <— > fk we get 

m . / 1 exp („„ _ J dh . PW) j |rr exp (w + ± _ ^±M)) . (16) 

We are now in the position to continue this result back to real values of the temperature by simply setting k = i[3. 
In this way we find the following self-consistent equation for the replica symmetric field distribution P{h) of a Levy 
spin-glass at inverse temperature f3 

m .J^ M -J M mJ^\rr^-^^^)) ■ (IT) 

V. SPIN GLASS TRANSITION 

From (fT7|) we infer that the paramagnetic field distribution, P(h) = 5(h), is always a solution. To test its stability 
we plug into the r.h.s. of (TTT| a distribution Po(h) with a small second moment, eo := / dhPo(h) h 2 <C 1, calculate the 
l.h.s. (to be denoted by P\(h)) by linearizing in eo and compare the new second moment, ei := J dhP\(K) h 2 , with 
eo. We find ei > eo, i.e. instability of the paramagnetic state, if the temperature T is smaller than a critical value 
Tf. a determined by 

fdrdf,,„ i„z ,o T(a + 1) ,a+l . f dr , 9 , . 

(Tf, a ) a = -J —\r\ a e* rr t & nh 2 f = - { - - ' cos(— tt) J ^-^ tanh 2 f. (18) 

This result for the freezing temperature is essentially the same as the one obtained by CB using the cavity method 
[HI]. Our somewhat more detailed prefactor ensures that the limit a — > 2 correctly reproduces the value Tf K = \/2 
of the SK-model [lfj ■ The dependence of TV a on a is shown in fig. [T] 
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FIG. 1: Freezing temperature T CtCl of an infinite-range spin-glass with Levy-distributed couplings as function of the parameter 
a of the Levy-distribution defined in ([2]). For the scaling of the coupling strength with TV as chosen in ([1} there is a finite 
transition temperature for all values of a. In the limit a — > 2 the result for the SK-model is recovered. 



The peculiarities of the spin-glass transition in the present system are apparent from the similarit y b etween (|17p 
and analogous results for strongly diluted spin glasses and disordered spin systems on random graphs [ll|, [2(| HH . To 
make this analogy more explicit we rewrite (| 1T[) in a form that allows to perform the s-integration to obtain 

J d=0 J i=l ' i=0 ' y ' 

°° LiW r d , dr df \ / 1 d \ 

= J2^T Tl( dh i p ( h i)^^\ r i\ a ^ iH ) ^(^--j^tanh-^tanh^tanh^)) . (19) 

(2=0 ' i=l " i=0 

This form of the self-consistent equation is similar to those derived within the cavity approach for systems with locally 
tree-like topology [ill . |20| and may also form a suitable starting point for a numerical determination of P(h) using 
a population-dynamical algorithm [2l| . 



VI. DISCUSSION 



Infinite-range spin-glasses with Levy-distributed couplings are interesting examples of classical disordered systems. 
The broad variations in coupling strengths brought about by the power-law tails in the Levy-distribution violate 
the Lindeberg condition for the application of the central limit theorem and give rise to non-Gaussian cavity field 
distributions with diverging moments. We have shown that it is nevertheless possible to derive the replica symmetric 
properties of the system in a compact way by using the replica method as developed for the treatment of strongly 
diluted spin glasses and optimization problems [2(| which focuses from the start on the complete distribution of fields 
rather than on its moments. 

Due to the long tails in the distribution of coupling strengths Levy spin-glasses interpolate between systems with 
many, i.e. O(N), weak couplings per spin as the Sherrington-Kirkpatrick model and systems with few, i.e. 0(1), 
strong couplings per spin as the Viana-Bray model. The majority of the N — 1 random interactions coupled to each 
spin are very weak (of order iV _1 / Q ). These weak couplings will influence only the very low temperature behaviour 
which may be expected to be similar to that of the SK-model. On the other hand the largest of N random numbers 
drawn independently from the distribution ([2]) is of order TV 1 /" [23| and hence every spin also shares a fraction of 
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strong bonds, Jy = Oil), which are for \Jij\ > 1/(3 practically frozen. With decreasing temperature a growing 
backbone of frozen bonds builds up that percolates at the transition temperature Tt a (l8| . The mechanism for the 
freezing transition is hence rather different from that operating in the Sherrington-Kirkpatrick model and resembles 
the one taking place in disordered spin systems on random graphs with local tree-structure. 
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